PERSISTENT CURRENTS IN INTERACTING ELECTRONIC 

SYSTEMS 



Georges Bouzerar and Didier Poilblanc 
Groupe de Physique Theorique, Laboratoire de Physique Quantique, 
Universite Paul Sahatier, 31062 Toulouse, France 

Persistent currents in disordered mesoscopic rings threaded by a magnetic flux are calculated 

\^ • using exact diagonalization methods in the one-dimensional (ID) case and self-consistent Hartree- 

0^ Fock treatments for two dimensional (2D) systems. For multichannel systems, a comparative study 

[ between models of spinless or spinfull (Hubbard) fermions has been done. First, it is shown that 

a purely one-dimensional model can not reproduce the correct order of magnitude of the observed 

currents. For 2D systems, going beyond first order pertubative calculations, we show that the second 

^ harmonic of the current is strongly suppressed in the case of spinless fermion models but significantly 

y—{ enhanced for the Hubbard model. This reduction (resp. increase) of the second harmonic is related 

^ \ to a strong increase (resp. reduction) of the spacial charge density fluctuations. Our work underlines 

' the important role of the spin degrees of freedom in the persistent currents. 
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I. INTRODUCTION 



The observations of mesoscopic currents in very pure metallic nano-structures was done in 
pioneering experiments In the first case, the experiment deah with the average current of 

^ a system of 10^ disconnected rings in the diffusive regime while, in the second, a single ring was 
G \ used. Although the existence of such persistent currents in small metallic rings was predicted 



^ • long ago [§,0, the magnitude of the observed currents is still a real challenge to theorists. Since 
^ ■ studies neglecting electron-electron interaction do not reproduce the order of magnitude of the 
O currents l^j-Q , there is a general belief that the interaction plays a crucial role in enhancing the 
^ current. It has been suggested that the currents could be enhanced to values close to the non 
disordered system case [pi?|,pT| . However, the exact role of the interaction in disordered systems 



is still unclear since treating interaction and disorder on equal footings is a difficult task. 

. _ I 

II. ID-SYSTEMS 

Let us first consider the simplest case of pure ID systems of spinless fermions. The hamil- 
tonian of the model reads, 

n = -t/2 {exp{i27r<^/L)c\cj + h.c.) +V nin^ + Y, wini (1) 

<ij> <ij> i 

where L is the size of the system and $ = where is the fiux through the ring and (pQ = —. 

The third term describes the disorder, where rii = CjCi and Wi are on site potentiel chosen 
randomly in [—W/2, W/2]. The nearest neighbour interaction has an amplitude V. 
If -£'($) is the ground-state energy the persistent current is given by, 

/ = ---. (2) 

An alternative to quantitatively estimate the persistent current is to calculate the Drude weight 
D, 



D 



L d^E 



(3) 



$m corresponds to the minimum of Using exact diagonalizations (ED) of small clusters by 

the Lanczos method [0, we have shown that, for strictly ID systems of spinless fermions, 
the effect of a repulsive interaction is to increase further the localization of the electrons and 
hence to decrease the value of the current. In fig.|l| we have plotted the Drude weight as 
a function of the inverse system size for two different band-fillings. At half filling and for 
W = 0, the Mott transition can be easily seen from this picture. The transition from metallic 
to insulator occurs at Vc = t (2t is the bandwidth). Away from half-filling and for W = 0, 
the system is always metallic and the Drude weight depends weakly on V. When disorder is 
switched on, for increasing repulsive interaction, we observe that D is reduced and, hence, the 
localization of the particles is increased. For V > t, the localization is both due to the Anderson 
localization and Umklapp processes generated by V. Note that, away from half-filling, we 
observe that the effect of V is very weak since Umklapp processes exist only in higher orders of 
the interaction. These results suggest that, for ID systems of spinless fermions, the interaction 
do not counteract the effect of the disorder but, on the contrary, increases localization. Using 



a Hartree-Fock approach, Kato et al. fl^ have obtained qualitatively good agreement with our 
ED results. 
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FIG. 1. Scaling of the Drude Weight at half-filling (right) and quarter-filling (left) vs inverse system 
size. The values of V and W (in units of t) are indicated on the picture. Each dot represents an average 
value over at least 300 configurations of the disorder. 



The case of attractive interaction is generically very different. Indeed, for some values 
of the interaction parameters, the persistent currents can be increased, even in the presence 
of disorder. For the sake of simplicity, we shall concentrate here on the half-field case. At 



fixed V, the system exhibits a Kosterlitz-Thouless-like transition |T5[ at a critical value of the 
disorder parameter pyc(^)- In fig-H we have plotted the phase diagram, calculated within a 
renormalization group (RG) approach, in agreement with ref. |[T^. This picture suggests that 



for —1 < y < —0.5 the system is metallic for sufficiently weak disorder. Note that W^""^ ~ 2.3t. 



For Y < — 1 and any filling, the system phase separates but this is of no interest for us in the 
present study. 
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FIG. 2. Phase diagram Wc(V') vs V obtained from an RG approach. 



We have compared the ED and the RG results ||T7| , p!8| and found a quaUtatively good 



agreement. In fig.y the Drude weight is plotted as a function of the inverse system size for 
various parameters V and W. For weak disorder W < t the agreement is clear. However, our 
numerical data suggest that the delocalized region is restricted to a smaller part of the phase 
diagram. This minor discrepancy might be due to the fact that the RG approach neglects 
Umklapp processes and is perturbative in V. 
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FIG. 3. ED results of the Drude weight vs 1/N for V < 0. Each dot represents an average value 
over at least 300 configurations of the disorder. 

The difference between the repulsive and attractive cases can be understood in the following 
way: in the case of repulsive interaction the 2fcj-fluctuations of the charge density are dominant. 
Therefore, the pinning by the impurity potential is strong and leads to a decrease of the 
localization length and, hence, to the reduction of the persistent currents. On the other hand, 
in the attractive case, the superconducting fluctuations are dominant. 

For spinfuU electrons (eg Hubbard model) Giamarchi et al. have shown that the in- 



teraction, on the contrary, enhances the persistent current. In this case, the increase of the 
current is closely related to the decrease of the spacial charge fluctuations or, equivalently, to 
the smoothing out in space of the charge density. Indeed, for the repulsive Hubbard model, the 
spin density wave (SDW) fluctuations are dominant in the ground state (GS). On the other 
hand, for the attractive case, the spin fluctuations become gapped and thus pinning is favored. 
However, it should be noticed that, when S> t, the charge degrees of freedom can be mapped 



onto a spinless fermion model with kf 2kf so that the 4kf component of the disorder is 



expected to become important in this hmit. At fixed W, we then expect a peak in the curve 
of D vs U. ED data on small Hubbard rings have indeed shown such a behavior. More re- 
cently, RG calculations taking into account Akf scattering have also confirmed this assumption 
pT| . Second order perturbation calculations are also in good agreement with the ED results 

Even though repulsive Hubbard interactions increase the persistent currents, 4/c/ scattering 
prevents the observation, in a strictly ID system, of amplitudes of the order of the experimental 
data. It should also be stressed that models of spinless fermions in the continuum lead to 
opposite results W^. Indeed, it has been argued that, in the case of ID continuous models, a 



repulsive interaction can increase the persistent current up to values close to the 'clean' case. 
These results emphasize the important role of both the spin and the nature of the interaction 
on the lattice. 



III. 2D-SYSTEMS 

As stressed in the introduction, real experimental systems are 3D anisotropic materials. 
Moreover, a diffusive regime can only be realized in d > 1. A complete understanding of the 
persistent current phenomenon then requires the study of multi-channel systems. In section HI, 
we shall present recent results concerning the effect of the electronic interaction for both spinless 
fermion and Hubbard multi-channel models [p3| , p^ . Note that we shall focus only on the second 
harmonic Ih2 of the current since (i) experimentally, in the multi-ring experiments, the current 
was found to be | periodic and (ii) theoretically, it is well understood that the ensemble average 



(over filling or disorder) suppresses the first harmonic of the current p5tP6| , P3 |. 

Diagrammatic first order calculations (the spin is here irrelevant) have shown that the 
persistent currents are increased by the interactions ||2^. More recently, Ramin et al. have 



numerically shown by a Hartree-Fock (HP) approach, that the first order correction of the 
persistent current was in agreement with the analytical treatment |^|. In both spinless or 



spinfull models, the second harmonic is enhanced. However, a nearest neighbour interaction 
between spinless fermions tends to decrease the value of the typical current while a repulsive 
extended Hubbard interaction enhances it. 

We have developped a self-consistent Hartree-Fock (SHF) method on finite lattices to treat 



the electron-electron interaction |23]. A systematic comparison with ED results on small lattices 



has revealed that this method was accurate even up to intermediate interaction strengths 
While, in the small coupling limit, the SHF method reduces to the HF method, at larger 
coupling it takes into account higher order powers of the interaction. Since the calculation 
is performed on (large) finite lattices quantum interferences due to the disorder potential are 
somehow treated exactly. It is important to note that this method is different from the usual 
perturbative scheme where corrections to the current due to the electronic interaction 
are calculated perturbatively. Through the self-consistency relation, our procedure includes a 
resumation of higher order terms which becomes essential at moderate interaction strength. 
However, the direct relationship with a standard perturbative expansion remains unclear yet. 
Our study could be applied both to single or multi-ring experiments. iFiom a theoretical point 
of vue, the difference relies simply in the absence or presence of particle number fluctuations. 




FIG. 4. {Ih2)dis ^ ^ function of V calculated within the SHF method for spinless fermions on a 
8x8 cylinder. An average over 1000 disorder configurations has been done. Comparison between first 
order HF (open symbols) and SHF results (full symbols) at half-filling. Circles correspond to W = 2>t 
and squares to W = At. 

Section III is organized as follows: first, we compare, for both spinless and Hubbard models, 
the first order correction (in the interaction) and the SHF correction to the second harmonic 
of the persistent current. As expected, the two methods are in very good agreement in the 
small interaction limit. For larger interaction strengths, higher order corrections (taken care of 
in the SHF method) become important for the case of the spinless fermion model while they 
remain small in the Hubbard case. Secondly, we show that, as the interaction is switched on, 
the decrease (resp. increase) of the persistent current is always closely related to the increase 
[lesp. decrease) of the charge density fluctuations. Lastly, a finite size analysis of the data of 
the second harmonic is performed. 

The starting hamiltonian reads: 



H = TLk + 'Hint + 'H 



'-des- 



(4) 



TCk is the usual kinetic part containing the flux dependance, TCdes is the term due to the 
disorder, 



(5) 



The sites i are now chosen on a L x L lattice with periodic boundary conditions in one direction. 
The disorder parameter W is restricted to a finite interval for which the (non-interacting) 



system is as close as possible to a diffusive regime [g^. Hint is the interacting part, where, in 
the Hubbard case, rii = + ni^. In the spinless fermion case. 



K 



int 



n 



(6) 



where i, j stand for nearest neighbour sites and V is the strength of the screened interaction. 
In the Hubbard case, the interaction part is local in space. 



'Hint = Uj2^i^ "-U- 



(7) 



In the SHF approximation, a mean-field type of decoupling is performed. 



which renormahzes the on-site disorder term (Hartree term) 5wi = Z^j^^i ^j('^j) ^he hopping 
term (Fock term) 6tij = Vij{cjc{). The quantities (rij) and (cj q) have to be determined self- 
consistently. Similarly, TC^^ becomes, 

nZt = t^E(K)^u + K)^iT - MM). (9) 

i 

For simplicity the calculations have been done in the paramagnetic sector (i.e. {uii) = (ni-^)). 
In this case the spin | and | are decoupled, 

'^fnt = ^'^int (10) 

(T 

where, 

ni^t = uY.iMni.-\M^). (11) 

i 

Note that, due to the disorder potential, the various mean-field quantities are space dependent 
so that a numerical investigation on a finite lattice is necessary. Similar equations also hold 
in the HF approximation but, in this case, the expectation values are simply taken in the 
non- interacting GS. 




FIG. 5. (a) (a) Ratio {Ih2) disi^^'^) / {^h2) disi^^'^) ^ ^ function of U calculated in the Hubbard 
model within both HF (open symbols) and SHF (full symbols) methods on a 8 x 8 cylinder. Circles 
correspond to VF = 3t at quarter-filling, squares to IF = 3t at half-filling and triangles to IF = 4t at 
half-filling, (b) Same as (a) at fixed W = 2>t and at half-filling for several system sizes. 



Since the current is periodic of period 1 (in units of $o), it can be expanded as a Fourier 
series, 

/($) =^4„s^n(27rn$), (12) 

n 

where Ihn are the harmonics of the current. In fig.§ we have plotted (//12) (() means average over 
disorder) at half-filling, for a 8 x 8 system, as a function of V. {Ih2) is calculated by averaging 



over at least 1000 configurations of tfie disorder. As expected, we observe for weak values of 
the interaction a perfect agreement between the SHF and the HF calculations. However, as we 
increase V, the SHF results show a strong reduction of the current while the HF results predict 
an increase. We also observe that the region of agreement between HF and SHF methods is 
reduced as the disorder increases. This means that, as the strength of the disorder increases, 
higher order corrections in the interaction become significant and, thus, a first order calculation 
is not sufficient. We will see later on that this reduction is related, as in the ID case, to an 
increase of the charge density fluctuations. Using the SHF method, we have checked that this 
effect persists for larger system sizes and away from half-filling. In conclusion, the repulsive 
interaction are detrimental to the persistent currents in spinless fermion models with short 
range interactions. Note however that, in a different regime, for much stronger disorder, an 
enhancement of the persistent currents by the interactions has been found pO |. 

However, a different behavior is observed in the case of the Hubbard model as seen in fig.^ 

showing the relative effect of the interaction on the current, i.e. (-^''^)('^'^\ Fig.^(a) reveals 

that both SHF and HF calculations predict an increase of the current and that the first order 
calculation gives relatively good results. It is interesting to note that the first order calculations 
always give higher values of the currents. In Fig.|^(b), we have plotted the same quantity, for a 
fixed density, but for different system sizes. It clearly indicates that, as the size of the system 
increases, the effect of U becomes stronger. 
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FIG. 6. Distribution -P(p) in both spinless fermion (a) and Hubbard (b) models, for a 10 x 10 
cylinder, and different band fillings n. The values of V, W and U are indicated in the figures. An 
average over 30 configurations of the disorder has been performed. 



At this point, these first results suggest the importance of the nature of the interaction. 
Indeed, in the spinless case, the currents are strongly reduced by the interaction while, on the 
contrary, the Hubbard repulsive interaction enhances the currents. In addition, it has been 
shown that in the spinless fermion case, higher order terms become rapidly dominant even 
for relatively modest values of the interaction strength so that a simple HF approach is not 
sufficient. 

Let us now try to build up a physical picture of this phenomenon. The distribution of 
the local charge density n-^ turns out to be directly connected to the localization. Considering 
many configurations of disorder (labeled by some integer k), we assume that the related local 
site densities nj' (where the subscript k stands for the disorder configuration) are independent 
realizations of a statistical variable p. In fig.^ we have plotted the distribution of the local 
density p for spinless fermion (a) and Hubbard (b) models. We clearly observe that the shape 



of the distribution changes with the interaction. In the spinless fermion case, two peaks clearly 
appear. This suggests the existence of two categories of sites with rather large and small 
charge densities, respectively. On the contrary, in the Hubbard case (b), we observe that the 
interaction has now the opposite effect i.e. the distribution shrinks around the average value. 
In other words, the interaction tends to homogenize the charge density in space. 
Let us now turn to more qualitative results. For this purpose, let us define. 



5pk 



\ 



1 

-piY.i{n>t) - n) 
^ i=i 



(13) 



and, 



5p 



N, 



^ cdes 



cdes 1^=1 



(14) 



where L is the length of the system, N^des is the number of disorder configurations, and (ji^) 
are calculated self-consistently. We expect, for big systems, 5pk to become independent of the 
disorder configuration. Results for 5p vs U or V are plotted in fig.|^ for various fillings. For 
spinless fermions, it appears that 5p increases with V at any fillings. For example, for n = 0.4, 
the width increases by almost a factor 4 at y = 0.8. However, in the Hubbard case, the effect 
of U is just the opposite e.g. the reduction of Sp{U, W) is larger than 25% at half filling. Note 
that, in general, the magnitude of the effect (relative increase or decrease with the interaction) 
increases with the density, as commensurability is approached. 
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FIG. 7. Ratio 6p{VU, W)/6p{0, W), (VU is V or U) as a function of the interaction parameter (U 
or V) for different filling factors n. The calculations have been done on a 10 x 10 cylinder with 30 
configurations of the disorder. 



To finish, we now turn to dependence of the width of the distribution 5p on the disorder 
parameter W summarized in Fig.|[ As expected, for vanishing interaction strength, the dis- 
tribution of the local densities increases almost proportionnaly with W. For finite interaction, 
Fig.^] shows that the sign of the effect observed depends only on the nature of the interaction 
(U or V) and not on the disorder. The magnitude of the increase/decrease seems to depend 
roughly on the combination ^ or U^. 
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FIG. 8. 5p{VorU, W) as a function of the interaction parameter (V or U) for different disorder W/t. 
The calculations have been done on a 10 x 10 cylinder with 30 configurations of the disorder. 



IV. CONCLUSION 

In conclusion, we have presented data beyond first order perturbative calculations. In ID, 
exact diagonalisation data clearly show that the electronic interaction enhances localization. 
In 2D, a self-consistent mean-field treatment of the interaction is used for the first time in 
the case of disordered fermions on a lattice. The SHF method has proven itself to be a good 
tool to study the effect of the interaction for d > 1. It provides strong evidences that spinless 
and spinfull fermions behave differently. The above SHF results show that the increase of 
the persistent current is strongly connected to the spin degrees of freedom. Although the 
calculations presented here concern 2D systems, we expect that, for bigger systems and higher 
connectivity, the pinning with the impurities becomes weaker. However, quantitative data are 
still missing for realistic mesoscopic sizes and more work is certainly needed to know whether 
such models can account for the correct magnitude of the experimentally observed current. 
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COURANTS PERSISTENTS DANS LES SYSTEMES DE FERMIONS CORRELES 



Les courants persistents dans des anneaux mesoscopiques desordonnes sont etudies a I'aide de 
methodes de diagonalisations exactes a ID et de methodes Hartree-Fock auto-coherentes a 2D. Dans 
le cas unidimensionnel d'un modele de fermions sans spin, les interactions electroniques augmentent 
toujours la localisation. En revanche, dans le cas de sytemes bi-dimensionnels de fermions avec spin 
1/2, le courant permanent est nettement augmente par les correlations electroniques. II est montre 
que ce phenomene est intrinsequement lie a la diminution des inhomogeneites de densite de charge 
dans I'espace provoquee par les interactions electroniques. Cette etude demontre I'importance cruciale 
des degres de liberte de spin. 



